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Abstract 

The distributions of two-jet event shapes contain information on hadronization in QCD. Near the 
two-jet hmit, these distributions can be described by convolutions of nonperturbative event shape 
functions with the same distributions calculated in resummed perturbation theory. The shape 
functions, in turn, are determined by correlations of momentum flow operators with each other 
and with light-like Wilson lines, which describe the coupling of soft, wide-angle radiation to jets. 
We observe that leading power corrections to the mean values of event shapes are determined by 
the correlation of a single momentum flow operator with the relevant Wilson lines. This generalizes 
arguments for the universality of leading power corrections based on the low-scale behavior of the 
running coupling or resummation. We also show how a study of the angularity event shapes can 
provide information on correlations involving multiple momentum flow operators, giving a window 
to the system of QCD dynamics that underlies the variety of event shape functions. In deriving 
these results, we review, develop and compare factorization techniques in conventional perturbative 
QCD and soft-collinear effective theory (SCET). We give special emphasis to the elimination of 
double counting of momentum regions in these two formalisms. 
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I. INTRODUCTION 



In perturbative quantum chromodynamics (QCD), scattering amplitudes computed with 
massless partons suffer from infrared and coUinear singularities, which, however, cancel in 
suitably-defined infrared safe cross sections. Event shapes, e(A^), are numbers associated 
with final states, N, defined so that the semi-inclusive cross sections with initial state / at 
center-of-mass energy Q, 



N 



are infrared safe for leptonic initial states (annihilation), and are factorizable into parton 
distributions and infrared safe short- distance functions for hadronic initial states. Infrared 
safety requires, in general, that the event shape e(A^) be equal for states that differ only by 
a rearrangement of collinear partons or by the emission or absorption of zero-momentum 
partons. 

Event shapes have been the subject of extensive study, partly as tools to explore the 
evolution of QCD dynamics from weak to strong coupling (see [l[ for a review). This paper 
will compare event shapes as treated in full QCD using factorization theorems and in the 
QCD effective theory, soft-collinear effective theory (SCET). In the process, we will derive 
a number of new results, relying on our ability to express infrared dynamics in terms of 
specific matrix elements in the full and effective theories. 

Although infrared safe, and therefore calculable self-consistently in perturbation the- 
ory, event shape cross sections at all but the very highest energies are not usually well- 
approximated by perturbative results alone. The relationship between perturbative and 
nonperturbative contributions to fully inclusive cross sections {e{N) = 1) was formulated 
early in terms of QCD sum rules [2| based on fixed order in perturbation theory and the 
operator product expansion, and later in terms of its high-order behavior [sj]. Eventually, 
it was realized that a similar analysis can be applied to many semi-inclusive cross sections 
that are infrared safe. 

The leading nonperturbative contributions to the mean values of many event shapes enter 
at the level of 1/Q in leptonic annihilation experiments, with Q the total center-of-mass 
energy. This dependence was predicted from the high-order behavior of QCD perturbation 
theory at infrared scales, using analyses based on the running of the coupling 0, H, 0, 0- 
These analyses naturally led to the proposal that coefficients of the 1/Q power are universal 
up to calculable factors. 

Power corrections in event shapes can be measured by comparing fixed-order perturbative 
predictions to data from a range of Q, and 1/Q corrections are readily seen in such compar- 
isons [;!]. For the mean values of many event shapes, these corrections are indeed universal, 
in a sense that we will review below. This is a striking prediction for nonperturbative effects, 
based on an analysis of perturbation theory at infrared scales. 

Beyond mean values, the distributions of event shapes require, as described below, con- 
volutions of resummed perturbative cross sections with nonperturbative "event shape func- 
tions" j^, [qI, a different one for each event shape. The event shapes summarize the leading 
power corrections for all the moments of the event shape in question. 

The roles and determinations of event shape functions are in some ways analogous to 
those of parton distribution functions. Parton distributions can be defined in terms of 
matrix elements in QCD, and these matrix elements have universality properties based on 
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factorization. As we shall review below, event shape functions may also be defined as matrix 
elements in QCD or SCET, which describe correlations between energy flow in different 
directions in the presence of lightlike color sources ^ . The universality properties of event 
shape functions, however, are certainly less well understood. 

These issues can be addressed in the terminology of perturbative QCD, using techniques 
of factorization and resummation. At the present time, however, a growing body of work in 
jet physics relies on the language of soft-coUinear effective theory to study power corrections 
[id . 11 | and perturbative resummations l^, T^]- With this in mind, we will use event shapes 
in the two-jet limit of leptonic annihilation to compare treatments based on factorized cross 
sections in perturbative QCD js], H, 0] and alternatively in SCET 0, We hope that 
this dual description will be helpful to some readers. Of course, even for this limited set 
of observables, a full comparison for both perturbative resummation and nonperturbative 
power corrections would require a lengthy discussion. In this paper, we will concentrate on 
power corrections implied by these closely-related formalisms. 

Not surprisingly, we will find that the two approaches are equivalent at the level of leading 
power corrections and shape functions that are associated with soft gluon emission. In both 
cases, infrared dynamics will be described by matrix elements involving Wilson lines. In 
particular, we will see that the soft shape functions of the factorization theorem can be 
defined to correspond directly to the functions that describe ultrasoft gluon radiation in 
these observables. 

Once we have discussed the formalism relating event shapes to matrix elements of mo- 
mentum flow operators, we will derive a number of results of phenomeno logical interest. 
The universality relations between different 1/Q corrections proposed in Refs. jl, [s], [g], 01 for 
mean values of event shapes can be derived from the field-theoretic shape functions, without 
invoking the dominance of single soft gluon emission or related assumptions on the univer- 



sality of the infrared running coupling [IJ]. In particular, we will derive the universality of 
leading power corrections to the average values of a large class of event shapes entirely from 
factorization and the boost-invariance of products of the relevant light-like Wilson lines. 

To specify the full event shape function for an arbitrary infrared safe observable, it is 
in principle necessary to know all the energy flow correlators. Relations between event 
shape functions, however, have been conjectured for a particular class of event shapes, the 



angularities [15|, [Ig, llTJ, which include the thrust and jet broadening. A scaling relation for 
the angularity shape functions follows from the assumption of negligible correlations between 
jet hemispheres for soft radiation. We will use insight gained from comparing factorized 
QCD and SCET to derive explicit relations between energy correlations and violations of 
the scaling rule, providing a set of measurements that relate directly to the correlations. 

Before going further, it is important to emphasize that the strength of nonperturbative 
power corrections to any observable depends generically on the definition of perturbation 
theory chosen for that observable. In particular, in observable-specific perturbative schemes 
for the coupling based on the method of effective charges [18], as reviewed in [l9|, the 
coefficients of power corrections to average values tend to decrease markedly compared 
to perturbative expansions in MS definitions. This method incorporates measurements of 
the observable in question directly into the renormalization scheme. It therefore builds 
more information into perturbative expansions for these observables than is possible in 
conventional schemes. Certainly a better understanding of the relationships between process- 
specific and process-independent approaches to perturbation theory would be helpful. In 
the discussion of this paper, however, we will assume that the strong coupling, as, is defined 
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in a process-independent fashion. 

Our discussion begins in the following section with a description of the class of event 
shapes that we study, and specifically defines the angularity event shape functions. In Sec. 
Ill, we review and relate factorization formalisms for event shapes in full QCD and SCET, 
exhibit the matrix elements that determine the soft gluon dynamics and define the event 
shape functions. We explore this relationship further in Sec. IV, with a comparison of 
the elimination of double counting through zero-bin subtraction in SCET 2^ and eikonal 



subtractions in QCD. In Sec. V we introduce momentum flow operators and draw the con- 
sequences of boost invariance for their correlations in soft functions. We show that power 
corrections associated with jet functions are subleading for a large class of angularities. 

In Section VI we apply the formalism of Sec. V first to the average values of event shapes, 
demonstrating the universality properties of these mean values. We go on to treat the 
scaling properties of angularities beyond their mean values, and we discuss the information 
on momentum flow correlators that is implicit in possible violations of the scaling rule 



proposed in Refs. [la, [17 



II. TWO-JET EVENT SHAPES 



In this paper, we study event shapes that can be expressed in the form. 



Q 



(2) 



where the sum is over final state particles, and the transverse momenta and (pseudo-)ra- 
pidities rji = In cot (6*4/2) are measured with respect to the thrust axis^ Two examples of 
these are the C-parameter and angularities (which include the thrust), which are expressed 



in terms of the rapidities as |16l . |21 



3 P.' 



cosh r]i'' 



Tn = — 



Q 



(3) 



In the limit e — > 0, the cross sections for all of these event shapes are dominated by final states 
consisting of two perfectly-collimated jets. For this reason, we refer to them as two-jet event 
shapes, and the limit e — > as the two-jet limit. Power corrections to the distributions for 
these event shapes enter at the level of l/(eQ)", star ting at n = 1, in addition to corrections 




suppressed by additional powers of Q 

Near e = 0, all such event shapes generate double logarithms in perturbation theory. 



which in many cases can be resummed to next-to- leading logarithms 22, 2^, 24] and beyond 
15|. A quantitative description of these distributions, however, requires nonperturbative 
input, which can be summarized in an event shape function, SNP,e- The physical cross 
section is then given as a convolution of this shape function with the resummed perturbative 
function js], 0], 



da{e, Q) 
de 



dCSNPACQ,^) 



dapTje - C,Q,A) 
de 



(4) 



^ The thrust axis is the choice of axis that minimizes the quantity on the right-hand side of Tq in Eq. ([3]) 
for a — 0. 
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The mass A, which we may think of as a scale comparable to, but larger than, Aqcd, repre- 
sents the boundary between perturbative and nonperturbative contributions. Below, it will 
appear as the renormalization scale for certain matrix elements that define the nonpertur- 
bative function. 

The nonperturbative function SNP,e{CQy ^) for event shape e is independent of the overall 
momentum scale, Q. Estimating Sjyp^eiCQy iii Eq. (jlj) from the plentiful data at the Z 
pole, for example, allows predictions of cross sections at lower and higher energies. These 
predictions successfully describe e"'"e~ data for the thrust and jet mass distributions for data 
over a wide range of center-of-mass energies js], 0, 25 



Our discussion below is focused on values of e that describe two-jet events, which dominate 
final states in leptonic annihilation. Extensive data on these shapes have been recorded, 
which may in principle be mined for information on the process of hadronization in QCD. 
Extensions to multijet events are possible, following Refs. 26, 13]. Certainly, the reasoning 
below will require further development for these cases. 



III. FACTORIZATION IN THE TWO-JET LIMIT 

In this section, we review basic results on the factorization of cross sections in the two-jet 
limit. We begin with an outline of the factorization analysis in full QCD, in the notation 
of Ref. 15], followed by a discussion based on soft-coUinear effective theory [l^, 11]. Our 



discussion will apply to both perturbative and nonperturbative contributions. 

Most of the results of this section have been given previously elsewhere, but we believe 
that a side-by-side presentation may help to shed light on both formalisms. In particular, a 
comparison of the formalisms will suggest the importance in both cases of the elimination 
of double counting. This will be the subject of Sec. IV, where we make use of the SCET 
discussion of "zero-bin subtractions" in Ref. [2C 



A. QCD matrix elements 



In the two-jet limit, the differential cross section (distribution) for a two-jet event shape 
e factorizes into a convolution of functions that characterize the jets with a "soft" function 
that describes wide-angle gluon emission 151 ]. 



= (^o{Q) / den den des S{e - e„ - - es)Jn{Q, e„) Jn{Q, en) Snn{esQ) + 0{e^) , (5) 

where Uo carries the overall dimensions, and can be defined as the Born cross section to 
lowest order in a^. In perturbation theory this cross section behaves as the order 1/e times 
logarithms as e vanishes, with contributions at order e° from wide-angle three-jet events. 

In Ref. 151], the perturbative resummation of the angularities was studied, with the 
functions in ([5]) defined in terms of QCD matrix elements. We may start by defining a set 



of path-ordered exponentials or Wilson lines. 



Pexp 



J — oo 



(6) 



where f labels a color representation. The vector S,ci which defines the path of the ordering, 
was taken in Ref. 15| to be off the light-cone, at least to start. 
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Following Ref. 



15 



, the jet functions of ([5]) are defined in terms of matrix elements 



2 (27r)^ 

g2 ^vT 



7^ 



$(?^(0)g(0) 



N, 



g(0)$gH0) 



X 6{ej^ - e{NjJ) 6 (Q - u;{NjJ) 6\nj^ - h{NjJ) 



(7) 



where c = n,n labels the direction of the jet as above. The jet functions are constructed 
from the squared amplitudes for the quark (or other partonic) field to produce states Nj^ 
with total energy Q, with a momentum whose direction n is in a fixed direction hj^, and 
with a fixed contribution ej^ <^ 1 to the event shape in question. The amplitudes are 
rendered gauge invariant by multiplying the partonic fields by ordered exponentials ([6]) in 
the directions and in the quark representation. 

To define the soft function, we introduce "eikonal" cross sections, 



0" 



(oik) 



(/^,e) 



^ $:(o|$i^)t(o)$(.)t(o) 



X (iV,J$(^)(0)$l'^)(0) O) (5(e-e(iV,0) 



(8) 



in which final states A^^ik are produced by products of Wilson lines in directions and 
n^, which we may take to be opposite-moving. The eikonal cross section, Eq. ([8]) must be 
renormalized, with scale /x. It provides a good approximation for soft radiation that is not 
collinear to these vectors, but it also contains collinear-singular radiation parallel to the 
directions of the lines. The collinear enhancements are already taken into account in the 
jet function, and the soft function S in Eq. must be defined in such a way as not to 
double-count these regions. In fact, it is easier than it might seem to avoid double-counting. 
This is because we can apply the same factorization to the eikonal as to the full cross section, 
factoring the same soft function S from a set of "eikonal" jet functions, which can themselves 
be defined in terms of matrix elements as 



^(cik) 



X (iV(-)|#(0)$(^H0) O) 5 (e, - e(iV(-))) 



(9) 



where the roles of the quark fields are taken by recoilless, lightlike Wilson lines. 

We first observe that there is a certain ambiguity in the separation of jet and soft func- 
tions. The ambiguity is exhibited clearly by a Laplace transform, where large u is conjugate 
to small e. Under the Laplace transform, da/de in Eq. ([5]) factorizes into a simple product 
of jet and soft functions. 



dee 



da{Q) 
de 



MQ) Jn{Q,v)Jn{Q.^)~S{Q/^) 



(10) 



Notice that dependence on the upper limit of the integral over e is exponentially suppressed 
at large u. The product on the right in (fTOj) can be treated in different ways, depending on 
the task at hand. The potential sources of double counting are eliminated in this case by 
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~S[Q/j,) = ——^ ^1^^^ . (11) 



defining S{Q/v) as the transform-space ratio of the eikonal cross section to the product of 
eikonal jets, 

In Eq. fill I) , the eikonal cross-section and eikonal jet functions approximate well the wide- 
angle soft radiation in the full cross section and in the jet functions, respectively. Thus the 
inverse eikonal jet functions in Eq. (|TT1) cancel the contributions of the soft radiation in the 
partonic jet functions in Eq. ffTOj) . When expanded in the coupling, the denominators can 
be reinterpreted as sets of nested subtractions. 

For the study of power corrections [9|], it will be more useful to implement a slightly 
different organization, found by simply shifting the factors Jn^^\Q,v)J^''^\Q^v) from the 
soft factors to the jets: 

a(Q, U) = (To(g) Jn{Q^ y)Jn{Q. ^^W^iQ, ^) ■ (12) 

where 



7 M - ■^'■W-'^) Jin „\ - ■^••W-") (131 



The eikonal subtractions serve the same role as above in Eq. (fTTj) . but now all double 
counting is subtracted from the jet functions. We will show in Sec. IIVI that this method of 



subtraction is directly related to the "zero-bin subtraction" scheme [20|] in SCET. But first 



let us review the factorization of jet cross sections from the perspective of SCET. 



B. Factorization for event shapes in SCET 

The SCET analysis begins with the expression for the distribution of event shape e in 
QCD (the "full theory" from an effective theory point of view), 

^ Y,mnQ)\0)L,W2T^)H\Q-pM)5{e-e{N)), (14) 



de 2g2 



N 



where = qV^q is the production current in QCD. The vector is the leptonic part 
of the amplitude for e"'"e~ --» qq. SCET organizes QCD in an expansion in powers of a 
parameter A ~ a/Aqcd/Q [281 . 2^. The modes in the effective theory we will use are 



collinear quarks and gluons, ^n,n and An^n-, and ultrasoft (usoft) gluons, Aus- These modes 
are distinguished by the scaling of the light-cone components p = {n ■ p,n ■ p,p±) of their 
momenta, defined with respect to light-cone vectors n,n = (1,0_l,±1). Ultrasoft momenta 
scale as pus ~ Q(A^,A^,A^). Collinear modes can have momenta with one of two possible 
scalings: Pn ~ Q{\'^,l,X) or p„ ~ (^(A"^, 1, A^), and similarly for n-collinear momenta pn- 
The theory with collinear modes with the first scaling is called SCETi, and with the second, 
SCETn tSfli]- The typical transverse momenta of the collinear particles in the final state 
determines the correct choice of scalings. (For example, very narrow jets with p± ~ Aqcd 
must be treated in SCETn.) We begin by matching QCD onto SCETj, and consider SCETn 
when we wish to account for narrower jets. 
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At leading (zeroth) order in the expansion in A, the full theory current matches onto 
the SCETi operators (3l| . 

(15) 
(16) 



where uj,uj' denote label momenta, and the jet fields Xn,n are given by 
Here Wn is a Wilson line of coUinear gluons, 

fO 



W„(z) =Pexp 



ds n ■ Anins + z) 



(17) 



Label momenta are the "large" pieces of collinear momentum. The momentum of a collinear 
mode splits into this label piece and a residual piece. More precisely, Pn = Pn + k, where the 
label momentum Pn = n ■ + p^ contains only the 0{Q) piece oi n ■ Pn and the 0{Q\) 
piece of p^ , and the residual momentum k is of order Aqcd in all components. A collinear 
field with a label u creates and destroys only modes with small O(Aqcd) fiuctuations about 
the momentum uo. 

The matching between the full and effective theories is performed by matching matrix 
elements at a scale /z: 



(^^)qCd(/^) = C(w,^';^)(j^ ,)gcET(^), 



(18) 



where the labels uj^uj' are summed over. The coefficient C(u;,u;'; /i) is the Wilson coefficient 
in this matching. The combinations of fields and Wilson lines in these expressions bear a 
close resemblance to the jet functions of Eq. ([7]). 

We can remove the coupling between ultrasoft gluons and the collinear fields in the SCETi 



Lagrangian via the field redefinitions [32 



A' V 



y'^w'y 



(19) 



and similarly for the n-coUinear fields, using the "outgoing" |33j Wilson line of ultrasoft 
gluons. 



Yn{z) = Pexp 



ig I dsn ■ Aus{ns + z) 





(20) 



For instance, the term containing the collinear quark-usoft gluon interaction becomes: 

Inin ■ Dusin I'jn ■ d^, (21) 

where iDus = id + gA^s- With this separation, we have established an effective theory that 
captures the dynamics in the full theory for the two-jet event shape distributions in the two- 
jet limit. The identification of jet and soft quanta precisely matches the "leading regions" of 
the full theory, as identified by analysis of momentum-space integrals for arbitrary diagrams 



3J|. This analysis, of course, is pro cess- dependent, and the vahdity of this SCET holds up 
to the same corrections that apply to Eq. ([5]), and is improvable by adding more jets, for 
example. 

The redefined jet production current becomes 



(22) 



8 



where the primes on the redefined fields x' refer to the use of ^' and W in Eq. (fT6l) . This 
leads to a factorization of the differential cross-section in Eq. (fT^ . analogous to Eq. ([5]), 

— = |C(Q,-g;/i)|^ dej aj{ej;fi)S{e- ej;fi), (23) 



where the function aj contains the coUinear fields and states: 

1 L'' 
2Q2 3 



Mej;/i) = ^^ \iNjMx:,Qr''Xn,-Q\0)\'W{ej-eiNj^NjJl (24) 



where is the spin-averaged, squared leptonic amplitude, and where the soft function S 
contains the usoft fields: 

5(e;/x) = -^Tr V |(XjT[y;,r^t] \o) \^)S(e - e{X^)). (25) 

We have split up the final state into the jets Nj^ - and the usoft sector imagining 
them to live in clearly separated regions. The notation in Eq. f l25]) setting the labels 00,00' 
of the coUinear fields to ±Q means that the fields must create a quark and an antiquark jet 
with total label momenta n ■ pj^ = n ■ pj- = Q, and = p^_ = 0. 

Before we proceed, let us consider briefiy the scale dependence of the jet and soft func- 
tions. The matching from QCD to SCET is done at the scale fi = Q, minimizing large 
logarithms in the Wilson coefficient C{Q, —Q; /i). The coUinear matrix element is naturally 
evaluated at a scale /ic ~ QX, and the soft function at fig ~ QX^, to minimize logarithms in 
those functions. The running between these scales and fi = Q is achieved by renormalization 



group evolution. The calculations required to do this in SCET are described in [12, ll3j. In 
this paper, we focus not on the details of the perturbative matching and evolution, but 
rather on the properties of matrix elements at the scale /i ~ Aqcd that must be treated as 
nonperturbative quantities. 

For observables dominated by fairly wide jets with p± ~ QX ~ a/Aqcd<5) the coUinear jet 
function at fic = QX may be evaluated in perturbation theory, leaving only the soft function 
at fig = Aqcd as a nonperturbative function. For observables in which narrow jets with 
P± ~ Aqcd contribute heavily, the coUinear matrix elements are naturally evaluated at the 
scale fic = Aqcd- Between the scales QX and QX^, the theory SCETj at the higher scale must 
be matched onto SCETn at the lower scale. The form of the coUinear operator (1221) remains 
the same; only the momentum scaling of the coUinear fields changes. A coUinear matrix 
element in SCETn is a nonperturbative quantity. In Sec. El we compare the relative sizes 
of nonperturbative corrections to event shape distributions from the soft and jet functions. 
For the observables of most interest to us, the power corrections from the soft function will 
dominate. 

Returning now to the jet function in Eq. (!24l) . we can further separate the matrix elements 
involving the jet quanta into the form: 

^Aej;fi) = ^^Tr [TlJ2{^\Xn,Q\NjJ{NjJXn,Q\0) (26) 



2Q2 3 

^r'^ E ^^--Q 1^^-^) (^^'^l ^k-Q |0) ] Wi^J - <Jn) - eiJn)) 



Nj 



N 
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The individual squares of matrix elements are closely related to the jet functions in 
factorization-based treatments of semi- inclusive cross sections, Eq. ([7]) [15]. In particular, 
the operator content, including the presence of the Wilson lines Wn, is essentially equivalent. 
The choice of vector (^^ in dZD) is somewhat different. Another difference is that here the 
collinear fields are fully separated from usoft partons, so that the sums over explicit jet final 
states |A^j„^) contain no usoft lines at all, while correspondingly, the virtual states that 
enter the jet matrix elements have no soft lines. In the soft function as well, jet lines appear 
neither in the final state nor in the matrix elements. 

In SCET, this bookkeeping is built in by a "zero-bin subtraction" prescription for collinear 
fields 13]. For each collinear field in the SCET Lagrangian, there is a sum over all label 



momenta, which must always be nonzero, to avoid overlap with usoft modes. Thus, regions 
of collinear loop diagrams or phase space integrals where a collinear particle can become 
ultrasoft must always be subtracted off. 

Each particle in the final state may be placed in a bin determined by its label momentum, 
leaving it with a residual momentum inside the bin, while those with zero label fall into the 
zero bin, and are assigned to the usoft sector The couplings of usoft particles to the 
states Nj^ - have already been factored out and accounted for in the soft function S. 

A natural question at this point is whether the jet and soft functions so defined are 
individually infrared safe. The answer is yes, because infrared safety does not require Lorentz 
invariance, only the hermiticity of the interaction Hamiltonian, and a sufficiently smooth 
shape function j35|. In our case, the cancellation of soft and collinear singularities in leptonic 
annihilation processes can be carried out for fixed values of one of the components of light- 



cone momentum [3J]. We thus expect infrared safety for each of the SCET jet and soft 
functions. 



IV. ELIMINATING DOUBLE COUNTING 



In this section, we continue our comparison of the SCET treatment of event shapes to 
analyses based directly on factorization 0, 15], concentrating on the elimination of double 



counting. The observations below are relatively simple, but to our knowledge they have not 
yet been made in the literature in this context. 

Proofs of factorization in perturbative QCD typically make use of Ward identities and 
subtractions to avoid double counting of leading configurations while maintaining infrared 
safety [l^, 36|. Unlike their SCET analogs, the matrix elements in Sec. IIIII above have no 



restrictions on their momentum integrals. The jet functions thus include a considerable 
contribution from what in SCET are classified as usoft gluons.^ The soft function 5* was 
correspondingly defined to avoid double-counting of collinear gluons by the soft function, and 
wide-angle soft gluons by the jet functions. As noted above, the soft function is constructed 
from Wilson lines in a manner precisely analogous to the SCET soft function, Eq. (!25ll . 
From Eqs. (fT2|) and (fT3l) . with double-counted soft modes subtracted from the collinear jets, 
we see that the eikonal cross section of ([8]) is exactly the same as the SCET soft function 



^ They also include quanta collinear to the opposite-moving jet. In SCET, these quanta are simply not 
present in the jet lines. In the factorized jet functions of Eq. ([7]), although present they do not produce 
collinear singularities because the Wilson lines are defined with respect to vectors that are not 
lightlike. 
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(125|) . except that the former does not include an exphcit hmitation to usoft final states, X„. 

The subtraction of double-counted soft modes from the jet functions is precisely what 
is achieved in SCET by the zero-bin subtraction for coUinear fields (cf. Eqs. (74) and (75) 
in Ref. [20]). Although there is a restriction to ultrasoft states in (12^ . the phase space 
integrals implicit in the sum still cover the entire region of momentum space, so that the 
eikonal soft function in QCD and the SCET shape function 0251) are practically equivalent. It 
is not inconsistent to integrate up to infinite momenta in a usoft integral. One is first taking 
the limit A — > 0, z.e. Q — »• oo, then integrating up to infinity the usoft momenta, which 
remain always formally much smaller than Q.^ It is important only that the soft functions 
for the observables we calculate are dominated by those particles which have "truly" usoft 
momenta. Below, we will examine the validity of this assumption for the event shapes we 
consider. First, let us elaborate on the connection between double-counting subtraction 
procedures in QCD and SCET. 



A. Zero-bin subtractions in factorization 

The leading-order (in A) Lagrangian of SCET contains several parts (see, e.g., Ref. (32|). 
First, the purely ultrasoft Lagrangian for quarks and gluons is identical to that in QCD: 

= qnsip^qus - ^ Tr G^^G,, + , (27) 

where the ultrasoft field strength, G^" = ^[D^^^, D'^^, is given in terms of the ultrasoft covari- 
ant derivative, iD^^^ = id'^ + gA'^^. The usoft gauge fixing and ghost terms are represented 
by /^Is^'. Meanwhile, the Lagrangian for coUinear quarks is 



in 



■ i^us + gn ■ + (n + 94tW^W\f^ + g4i^^,) 



(28) 



where V and V± are label momentum operators which pick out the 0{Q) piece of the n ■ p 
component and the 0{QX) piece of the p± component of the coUinear momenta. Here W 
is the Wilson line defined in Eq. f ll7l) . organizing coUinear gluons n ■ A, each of which are 
leading power in A [33]. There is also an implicit sum over all labels and conservation of 
label momentum in each term [37]. Finally, the coUinear gluon Lagrangian is 



(29) 



where iV^ = + + ^in ■ D^g, and Cf^' contains the coUinear gauge fixing and ghost 
terms. 

The component n ■ A^^ of the usoft gluon field appears in the coUinear quark and gluon 
Lagrangians above. The field redefinition (fT9!) removes this interaction, through the identi- 
ties in ■ D^,Y^ = and W = Y^Wf^-^Yn, where iPj,) = + + ^in ■ d. This gives the 



'us-^n ~ " "'^^'-^ ''^ ~ n ''■^(0) 

redefined coUinear quark Lagrangian, 



We thank A. Manohar for discussions on this point. See related remarks in Ref. 11 1 
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ri (q) _ f ' 



2 



(30) 



and coUinear gluon Lagrangian, 

4^^^ = ^ Tr { + gA'l^, ^V^^,^ + gA'^,?} + C'f' • (31) 

It is in the theory described by the redefined Lagrangian that we derived factorization 
for the two-jet event shape distributions. As given in Eq. (12^ . the full distributions are 
convolutions of soft and separate jet functions for the two jets, which are easily abstracted 
from Eq. ( l26l) . For example, in terms of the coUinear SCET quark field, S,'n in (|30ll . we can 
define a jet function J^(e), as 

Jn{e) (I) = 5^ |(iV.J [«k |0)|^,/(e-e(iV,J), (32) 

where the total label momentum of the coUinear fields must equal the label momentum 
pn of the states A^j„, which are free of quanta in the usoft sector. The right-hand side is 
proportional to the matrix ^/2, as shown, and o/^(e) is a remaining scalar function. Here 
we include a subscript on the matrix elements to emphasize the Lagrangian in question. 
This jet function contains contributions only from "truly" coUinear particles with nonzero 
label momenta, and all usoft particle contributions are accounted for in the soft function 
5*. Restricting coUinear label momenta to nonzero values avoids double-counting of usoft 
contributions in these two functions. It may be convenient, however, to allow integrals over 
coUinear momenta to include the "zero-bin" . One must then define a prescription to subtract 
it out again. This can be achieved by the zero-bin subtraction procedures illustrated by 



Manohar and Stewart in Ref. [20], or by the method of factoring out "eikonal" jet functions 
described by Ref. [isl]. 

In order to establish the connection between the zero-bin and eikonal jet methods for 
subtraction, let us take another look at the coUinear quark Lagrangian in SCET, now with 
sums over labels explicit: 



C'}'\x) = J2 «^^''"'^-^Cp'(x)[m-9 + ^5^e-^^--n-<,^(x)]|^^^^ (33) 



We have dropped for now the terms with j^-^ gluons, since they are irrelevant for the remain- 
der of this discussion"^. The tildes denote label momenta containing only the large 0{Q) 
and 0{QX) pieces, = n ■ p^ + p^, and the remaining x dependence of coUinear fields 
fiuctuates on a scale described by residual momenta k of order QX^. The sums over coUinear 
labels are restricted to nonzero values to avoid double counting the usoft modes, and we have 
shown explicitly the phases that enforce label momentum conservation. As in the zero-bin 



prescription of Ref. [20[ , the sums over labels can be allowed to include the zero bin if we 
add a term that subtracts it out again. Let us worry about this only for the gluon field, 
since the interaction of usoft quarks with coUinear fields is subleading in A. These steps can 
be achieved beginning with a Lagrangian C" that includes zero-bin coUinear gluons: 



^ The perp gluons in the zero-bin overlap with usoft gluons, but these interact with coUinear fields only at 
subleading order in A. 
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We then perform a field redefinition on (nonzero-bin) collinear fields: 
^" —Tl'^f' A" — Tl'f A' TT W" — TT^W'TT 



in 



a + ^7$^e-^^"-n-<^(x)]|c,p(^ 



where Un is a "zero-bin" Wilson line: 



Un{x) = Pexp 



tg I ds n ■ A'^ ^{ns + x) 



(34) 



(35) 



(36) 



mimicking the field redefinition (fT9l) with the usoft Wilson line. The Lagrangian obtained 
from C"^"^^ by the redefinition (1351) is precisely C'J''^\ with the q = term subtracted off. The 
redefinition of the gluon field removes the zero-bin gluons from the corresponding collinear 
gluon Lagrangian. Note, however, that it does not yet remove zero-bin n ■ fields in W". 
Thus, the Wilson line W is not yet W, when the latter is defined to be free of zero-bin 
gluons. 

The jet function J^(e) of Eq. fl32l) is, as indicated, calculated in the theory described by 
C'^ in order to avoid double-counting the usoft modes. However, if we for convenience decide 
to include the zero-bin in all the collinear momentum integrals, we are actually calculating 
a slightly different jet function, J„, in the theory described by C'^. Via the field redefinition 
( 135|) we can move back to the Lagrangian C'^, but then we must include the Un Wilson line 
in the current: 



JJe) 



i 



Y.\^NjM'nW'X. |0)|^„5(e-e(iV,J) 



Y.VNj^\ &'X.Un |0) 5{e - e(iV^J) 



Nj 



(37) 



In the Lagrangian however, the zero-bin gluons do not interact with the collinear fields 
at all. Once the n ■ A^^ zero-bin gluons have been removed from W , this jet function 
factorizes into "purely" collinear and "zero-bin" parts. To accomplish this, let us also split 
the collinear Wilson line W^^ into its zero-bin and purely collinear parts. 



W'niz) 



where 



VLn{z) = Pexp 



ds n ■ Anfl{ns + z) 



(38) 



(39) 



and W'^ is the collinear Wilson line with restrictions to nonzero label momenta on all gluons. 
Then the alternative (£") jet function factorizes into: 
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where is the jet function calculated in the purely-coUinear theory described by and 
Jn^^^ is the eikonal jet function. The collinear fields in can only produce gluons with 
nonzero label momenta, while the Wilson lines [/„ in Jn^^^ produce gluons in the zero 
bin. It is closely related to the eikonal jet function defined above, in Eq. ([9]). 

By taking the Laplace transform of the jet functions in this convolution, we may solve 
for the purely collinear jet function: 

which is precisely Eq. (fT3l) . This relation tells us that we may calculate the jet function on 
the left, J^, by calculating instead the jet function on the right, Jn, which include the zero- 
bins in collinear integrals, if we then subtract out the double-counted usoft contributions 
with the eikonal jet function, which is precisely the prescription of Ref. [l^. We emphasize 
that the alternative SCET jet function, J„, is directly analogous to the pQCD jet function 
defined in Eq. ([7]). There remain technical differences associated with the states \NjJ in 
([7]), which can also contain energetic lines collinear to the opposite-moving jet, for example. 
As noted above, the effect of such lines is perturbatively calculable, so that these remaining 
differences are perturbatively calculable. We have thus shown its equivalence to the zero-bin 
subtraction procedure of Ref. [20| in SCET. 



B. The ultrasoft region in SCET 

For the purposes of our discussion below, we need to consider in some detail the precise 
definition of our usoft region. In each of the two-jet event shapes of Sec. 2, the contributions 
of particles in the far forward and backward regions are suppressed exponentially in the 
absolute value of their rapidity, and hence as a power of their larger light-cone momentum 
(order Q). The contributions of jet-like particles at low (order A) transverse momenta are 
suppressed by a power of Q compared to those emitted with similar low transverse momenta 
at wide angles. For this reason, to leading powers in l/{eQ), we need not include separate 
nonperturbative functions for the jets |5|, which may be expected to enter beginning at 
the level of l/{eQ^~^^), where b depends on the event shape at hand. For the angularities, 
b = 1 — a. We will give a more formal argument for this result in the next section. The soft 
shape function organizes all powers like l/{eQ)"', while neglecting all further suppression by 
powers of Q. 

Our analysis of power corrections from the soft function will depend crucially on the 
boost invariance of the usoft Wilson lines appearing in S, and, in the sum over usoft states 
in (!25l) . we integrate over all momenta. With this in mind, we will want our ultrasoft region 
formally to include all gluons with small transverse momenta but boosted to arbitrarily large 
rapidities. The light-cone components of their momenta are then, strictly speaking, larger 
than the typical usoft scaling QX^. Formally, the usoft modes in SCET do cover all momenta 
up to the scale fi = Q, although the effective theory Lagrangian is a good approximation to 
the dynamics only of the "truly" usoft particles. As long as we pick observables to which only 
these truly usoft particles contribute significantly, it is safe to include the highly-boosted 
particles as well in the usoft sector of the theory. Indeed, for event shapes such as the 
angularities, the contribution of particles in the far-forward and far-backward regions are 
p ower-suppressed . 
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Consider for example the shape function (l25ll for the angularities as defined in Eq. ([3]). 
Do the usoft modes in SCET correctly describe the dynamics of all the small transverse- 
momentum particles that make a non-negligible contribution to the observable r^? Let us 
say that if the exponential factor, exp[— \rii \ (1 — a)], for particle i is of order A or smaller, 
then its contribution is negligible. Then the largest-rapidity particle making a non-negligible 
contribution to the event shape has: 

Ir^il ~ --^InA, (42) 



or, defining the rapidity as ?7 = ^ ln(?T, ■ k/n ■ k), 



1 



( n ■ k n ■ k\ , ^ , 

max — r ~ ■ (43) 

\n ■ k n ■ k J 

Now, because n ■ k n ■ k k\, and for usoft particles, k± ~ QX^, Eq. ( H3|) implies that usoft 
particles with the larger of their light-cone momenta up to the order QX i-" contribute 
non-negligibly to the event shape r^. For a < 1/2, this light-cone momentum is still smaller 
than the corresponding component of a coUinear momentum, which is order Q [e.g. for the 
thrust, a = 0. With Q ~ 100 GeV and Aqcd ~ 1 GeV, we have A = 0.1, so the largest 
usoft light-cone momentum that contributes to tq is 10 GeV, still well below Q.) As long 
as this hierarchy of scales holds for the large light-cone components of usoft and collinear 
momenta, the only component of usoft gluons that interacts with collinear modes in the n 
direction in the leading-order SCET Lagrangian is the n ■ Aus component {n ■ A^s in the 
n direction), so that these interactions can be removed by the field redefinitions with the 
Wilson lines Yn^. This guarantees the form of our usoft shape function fl25l) . and ensures 
that only the "truly" usoft particles contribute to the sum. We can extend the range of 
allowed values of a beyond a < 1/2 (but only up to a < 1) by relaxing our criterion for 
the size of "non-negligible" terms {i.e. allowing exp[— |?7j| (1 — a)] to be larger than A but 
smaller than 1). For the C-parameter, we need 1/ coshr^.^ < A, which for the values chosen 
above translates to n ■ k,n ■ k < 20 GeV. 

Thus, in the following we may safely incorporate the power-suppressed contributions of 
the very far-forward and far-backward radiation in the SCET shape function, Eq. (l25ll . and 
identify it with the eikonal cross section, Eq. ([8]), evaluated at the corresponding scale. This 
means that in both SCET and full QCD we may treat the sum over states in the shape 
function as boost-invariant. It is this result that will lead us to demonstrate universality 
properties below. 



V. IDENTIFYING POWER CORRECTIONS 



We are now ready to identify the power corrections that arise naturally when the soft 
and jet functions are evaluated at scales of order Aqcd- In doing so, we set aside issues 
of perturbative resummation and of matching, treated in full QCD for the angularities in 
Ref. |15|, and very recently in SCET for the closely-related jet cross sections by Refs. H, 13 



In the discussion of this section, we will find useful a variant of the energy flow operators 
introduced in in Refs. j38| and applied in this context by Ref. This operator is clearly 
also closely related to the energy-energy correlations of Refs. |39|. 
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As mentioned earlier, the typical transverse momenta of the collinear particles in the 
jets which contribute to a given observable determine whether they should be treated in 
the theory SCETj, in which collinear momenta scale as Pc = {n ■ Pc, n ■ PcP't) ~ Qi^^y I5 ^) 
or (5(1, A^, A), recalling that A ~ v^Aqcd/Q- The typical virtuality of such particles being 
pI = QAqcd, jet functions in this theory can be calculated perturbatively. However, some 
event shapes may weight much narrower jets more heavily, in which jet constituents with 
transverse momenta of order QX^ ~ Aqcd become important. These degrees of freedom 
must be treated as collinear particles in SCETn, in which collinear momenta scale as pc ~ 
(^(A"^, 1, A^) or Q(1,A^,A^). These particles have virtualities of order pi ~ ^qcd^ ^o 
give rise to nonperturbative effects in addition to those from soft particles. For such event 
shapes, nonperturbative power corrections to the jet functions may compete with (or even 
dominate) those in the soft functions. 

We will now show how these results can be justified for the event shapes in question. 
We will give our arguments in terms of SCET matrix elements, keeping in mind that they 
can be presented in terms of matrix elements in full QCD in a similar manner. For the 
C-parameter and angularities Ta with a < 1, the dominant power corrections (of the order 
Aqcd/Q) will come only from the effect of usoft particles whose momenta are of (^(QA^). 
Power corrections from collinear particles will be found to scale as (Aqcd/Q)^ which 
then dominate for a > 1. However, for a > 1, there are also other power corrections, for 



example, due to the shift in the thrust axis itself caused by the soft radiation |15|, |23|, 140 
The inclusion of these effects, while necessary for a complete treatment of power corrections 
to Ta with a > 1, is outside the scope of this paper. 

Consider the distribution of an event shape of the form in Eq. ([2]), given in SCET by 
Eqs. fl23l - [25|) . The collinear cross-section fl2^ is, writing out the general event shape of 
Eq. ([2]) in the delta function explicitly, 

^j(ej;/ic) = ^ E \(NjnNjJXn,Qr^Xn,-Q\0)f{f^c)s{ej-^ Yl \P^\feind), 

(44) 

while the soft function is 

5(e;/x,) = i-Tr ^ |(X„| Yj'n |0)|' (/i.)5(e - ^ ^ |k,^| f^), (45) 

where we have now chosen to denote explicitly the dependence of the jet and soft functions 
on the scales /Xc/^s- Also, we have suppressed the factor associated with the leptonic part, 
and we have removed the time-ordering operator that was in the soft function in Eq. (!25|) 
by using the Wilson line Y^, where the bar denotes the anti-fundamental representation of 
SU{Nc) [iH. For event shapes such as Tq for a < 1, the collinear scale /ic can be chosen at 
a perturbative scale /Xc ~ QA, and we are in SCETj. For a > 1, the event shapes pick out 
narrower jets so that the collinear scale is determined to be of order /ic ~ Aqcd, putting us 
in SCETn, where the jet function is nonperturbative. 

We may express the delta functions in Eqs. (jH]) and (1451) in operator form by making 
use of a transverse energy flow operator, defined by its action on states N: 

St{v) \N{h)) = K\ ^iv - V^) \N{h)) , (46) 
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where the sum is over the particles i in state A^. This is equivalent to the energy flow 
operators discussed in Refs. [o], 38, 31 • In terms of this operator, the collinear and soft 



functions fl44ll45p can be written as 

1 - / 1 1"°° 

^j(ej;/ic) = ^ J2 {^\Xn,-QT^Xn,QS(ej-- dr] fMSAv)) \Nj„Nj^) 



X {Nj„Nj,,\xn,QT^Xn,-Q\0), (47) 



and 



^e(e; /X.) = ]^ Tr J] (0| yIy^6 " ^ feiv)£Av?j \Xu) (XJ Y^Yn |0) . (48) 

We can expand the delta functions in power series to identify the power corrections. If we 
flrst factor out the overall, canonical factor of 1/e, shared with perturbation theory, and 
assume that the matrix elements are of the order of the momentum components of the usoft 
gluons, QX^ ~ Aqcd, we derive a power series in AQCD/(e(5). Indeed, the purpose of event 
shape functions is to organize all terms in this series when Aqcd/IcQ) ~ 1 and all such power 
corrections are comparable. These power corrections are particularly clearly exhibited by 
Laplace transforms, Eq. (ITU]) , of the soft function at low scales jof. 



Se{i^]fis) = / deexp[-ue] Se{e; fis 



-^Tr(0|4F„texp 



g j dr] fe{r])ST{v) 



Y^Yn |0) , (49) 



where we have summed over the complete set of intermediate states in (1481) . as argued in 
Sec. IIVBI Expanding the exponential, we flnd a series in powers of the Laplace variable v 
divided by Q. 

For the soft function, according the discussion in the previous section, the sum over usoft 
states is unrestricted, as is the integral over rapidities inside the delta function. In the 
collinear function, choosing ej to be close to the two-jet limit ej = 0, or specifying a jet 
deflnition to pick out two-jet events, restricts the phase space integrals in the collinear cross- 
section to those with large rapidities, effectively limiting the range of the rapidity integral 
as well. The rapidities rj can be written in terms of the light-cone momenta of flnal state 
partons, f] = ^ ln{n-p/n-p). For a usoft parton, the ratio n-pus/n-Pus ~ 1? as all momentum 
components are 0{QX'^), so |?7| ~ 0, while for collinear partons in SCETn, one light-cone 
component is 0{Q) while the other is 0{QX^). Thus, n-pc/n-pc ~ or A~^, so r/ ~ ± In A^. 
Consider what this implies for the collinear and soft functions in the case of the angularities. 
The function /e(e) for e = Ta is fraiTa) = e"''''*^^"'^^. In the usoft function, this factor is of 
0{1). In the collinear integral, the phase space restrictions limit the rapidity integral to 
\ri\ > ln(l/A^), so that the collinear function is effectively 

^ Nj^Nj- ^ "^I'll? 

X {Nj„NjJxn,Qr^Xn,^Q\0). (50) 
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Although we cannot compute these nonperturbative matrix elements at the scale /ic ~ Q^^, 
we can estimate their dependence on A from dimensional analysis. Matrix elements of powers 
of the operator Sxir]) in coUinear states in SCETn should vary as corresponding powers of 
QX^. Similarly, each rapidity integral should behave as X^^^~"'\ Combined with the factor 
1/Q in front of the rapidity integral, power corrections to the coUinear jet function occur 
as powers of X^~^"'/Ta = (l/ra)(AQCD/Q)^~"- Correspondingly, in Laplace moment space, 
this becomes a power series in z/(Aqcd/Q)^ The latter is also the only argument for the 
jet function that serves as a boundary condition in the perturbative QCD resummation of 
Ref. [3].^ As long as a < 1, we may consider these to be sub leading compared to the power 
corrections of the soft function, which are powers of Aqcd/Q- For a > 1, we must take them 
into account, along with the recoil corrections mentioned above [l^, 17. 4ll . 

From now on, we consider only observables that pick out jets with typical transverse 
momenta well above the nonperturbative scale. In the language of SCET, this allows us to 
work in the theory SCETj and consider power corrections only from the soft function. 



VI. MOMENTUM FLOW OPERATORS, UNIVERSALITY AND SCALING 
A. Nonperturbative Universality from Perturbative QCD 



A striking prediction from the analysis of event shapes in perturbation theory, including 
those given in Eg. (3|), is the universality of power corrections to their mean values [sl, [g], 0, 



i, E, m M, MM, M 



(e) = (e)pT + c, 



A 

Q 



(51) 



In this expression, A a universal parameter and Cg is a calculable coefficient that depends on 
the observable, as we shall see below. The same reasoning that leads to fl5T|) . when applied to 
the event shape distributions, produces a shift in the resummed perturbative cross section. 



da 
de 



PT 



NP 



da 
de 



A 



(52) 



PT 



These relations were derived in Refs. [6|, |7|, |42| from the assumption of a "dispersive" repre- 
sentation for a.sifi'^) considered as an analytic function of the scale fi, and in Refs. jl| they 
were abstracted directly from the form of resummed perturbation theory. 

A more general approach jo], [l^, 17] replaces the shift of Eq. fl52l) by a convolution with 
a shape function defined as above, which reduces to a product in Laplace moment space, 
Eq. flTU]) . As we have noted, these shape functions are all different, but for the angularities 
a generalization of the universality of Eq. fl^T]) has been suggested, in the form of a scaling 
relation. The Laplace-transformed shape function for angularity distributions arising from 
resummed perturbation theory at next-to-leading logarithm (NLL) 
scaling with the parameter a: 



17| displays a simple 



1 



1-a 



n=l 



Q 



(53) 



See, for example, Eqs. (67) and (74) of [15 1 
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where A„ is independent of a. If we keep only the hnear, v /Q, term in the shape function, its 
inverse Laplace transform gives a delta function, which in the convolution of Eq. ([5]) leads 
immediately to the shift of Eq. (!52l) . As noted above, we limit our attention to angularities 
for a < 1. The values of the coefficients of Eq. of course, must be abstracted from 
a combination of experiment and resummed perturbation theory. 

Event shape functions derived from resummation organize all corrections in zz/Q that are 
implied by perturbation theory. Formally, the coefficients A„ are given in the NLL resummed 
cross section by 



z/V 1 dpi 
nn\ \ Q) J pI 



A„ = r^(-7^) I ^P^rAia^ipT)) , (54) 





where A{as{pT)) = Ci(tts/vr) + ... is the residue of the 1/(1 — x) pole in the splitting 
function for the parton, i = quark or gluon, that initiates the jet, and k, is an infrared 
factorization scale. In this picture of power corrections, the coefficients A„ are independent 
of a. The coefficient of the lowest power, n = 1, is equivalent to an integral over the running 
coupling, defined in a scheme where the coupling incorporates all higherpowers of A{as) in 



MS j45| . This approach generalizes the dispersive treatment of Refs. [7|, |42, |43| to higher 
power corrections, but shares with it a reliance on (exponentiated) low-order gluon emission. 

Analyses based on a dispersive coupling or on resummation rely to a greater or lesser 
extent on the kinematics of single soft gluon emission in the final state, and the universality 
relations follow from the boost invariance of these emission cross sections. The "Milan 



factor" [43| of the dispersive approach accounts for effects at next order in as, where boost 
invariance and the resultant universalities can be maintained. We are about to show that 
the boost invariance of Wilson lines in the soft shape function Eq. fl2^ is by itself enough to 
prove the universality relation for the mean values without further assumptions. In Refs. 



and [4^ the role of energy flow was explored in a manner closely related to our discussion 
below. 



B. Universality in Average Event Shapes from the Soft Function 



We continue to limit our attention to event shapes for which the dominant power correc- 
tions come from the soft function. For the shape function in the form given in Eq. (jUD, the 
operators in the matrix element no longer contain any reference to the flnal state X^, so, as 
in Eq. P9|) . we may perform the sum over intermediate ultrasoft states, leaving 



^e(e) = -^Tr 



' yiy^i^ ( ^ - Q / dv^T{v)fe{v) ) ynyn io) 



(55) 



From now on we drop the explicit dependence of the soft function on the scale /i^. In (l55l) . 
we insert factors of U{A{7]')yU{A{r]')) = 1, implementing a Lorentz boost of each operator 
in the 2;-direction with a rapidity rj'. The vacuum |0) is invariant under Lorentz boosts, and 
the Wilson lines are also invariant: 



U{Aiv'))yn{0)U{A{r]')y = f/(A(r/'))Pexp 



ig dsn ■ A^sins) 



U{A{r^')V (56) 



Pexp 



ig I ds an ■ Aus{ans) 





yn{o), 



19 



where a = e ^' , as n ^ an and n ^ a ^n. (This is also known in SCET as type-Ill 
reparametrization invariance [i^.) The only change is in the operator Sxif])'- 

U{A{v'))ST{v)U{A{7i')y = St{v + V') , (57) 

which follows from the defining relation for the St operators, Eq. f j46l) . Thus, the argument of 
the operator £'t('7) in the shape function in Eq. ( l55l) may be shifted to any value of rapidity, 
St{v) ~^ ^t{v + v')- At this stage, this does not yet allow us to perform the rapidity 
integral of fe{v) inside the delta function. Thus we do not find that the leading power 
correction simply shifts the argument of the perturbative event shape distributions, as the 
delta function is a highly nonlinear function of the energy flow operator and sits sandwiched 
between Wilson lines in the matrix element. If we do neglect correlations between these 
operators, we derive a delta function for the shape function, and reproduce the shift in the 
distribution, Eq. ([52]) 

The boost property fl57j) of a single operator, however, gives a strong result when applied 
to the first moment of an event shape distribution . Taylor expanding the delta function 



in Eq. (I55p (which is valid if we integrate the distribution over a sufficiently large region 
near the endpoint), we find 

^,(e) = (5(e) - 5'(e)i J dv fe{v)^ Tr (0| FtF„t£:^(^ + |0) + ■ • • . (58) 

Recalling the boost properties of the Wilson lines and the energy fiow operators St{i]), we 
are free to choose any value for rj' in this expression. Then, choosing rj' = —rj, we find that, 
remarkably, we may take the matrix element of the St operator out of the integral over i], 
leaving the result 

S,{e)=S{e)-S'{e)c,^ + -- - , (59) 
where the coefficient Ce is given by the integral, 

POO 

Ce= / dr]Mri), (60) 



oo 



and the universal quantity A is 

^ = 1 Tr (0| YiY^STrnJ"^ |0) . (61) 
For the C-parameter and angularities r^, the integrals of the corresponding weight functions, 

over all rapidities give the coefficients, 

2 . , 

cc = 37r, Cr, = J— ^. (63) 

When convoluted with the perturbative distribution, 5'e(e) reproduces the universality re- 
lations of Eq. (|5T|) for the first moments of the distributions. We have thus established 
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these results without appeahng to a one-gluon or related approximation. All higher-order 
corrections due to multiple-gluon emission separate from the observable-dependent factor 
Ce, which can be computed in a "naive" fashion [1] as in Eq. (l60l) above. 

The result for the C-parameter may be extended to a larger class of related event shapes 
by defining functions, fc^ [f]) =3/ cosh" rj, by analogy to the angularities. The integral over 
rapidities of this function gives the coefficient cc^ = 3B{a/2, 1/2), where B{x,y) is the beta 
function. In like manner, various new event shapes may be defined by appropriate choices 
for the function fe{v)- 



C. Angularity Distributions and Momentum Flow 



The expression ( !55|) for the shape function in terms of energy fiow operators enables us 
to put the power expansion of Eq. (153|1 into a more general field-theoretic context, and to 
discuss the possible significance of scale breaking. 

Let us compare Eq. (l53ll . derived from resummed perturbation theory, with the Laplace 
transform of the corresponding shape function in Eq. (!55ll (i^ ]. This is given by 



Tr(0|4F„texp 



c 



Q 



YnYn |0) , 



(64) 



which can be re-expressed as an expansion in cumulants. 



In 



Sa{u) 
^„(0) 



oo / 

y-i-- 



n=l 
oo 



^,,,e-lil<i-»)fy(,,) 



n=l ^ ^ ' 



(65) 



Here, and below, in the cumulants the Wilson lines Yn and Yn are understood. With this 
normalization, the coefficient Ai{a) for the angularities is related to the universal coefficient 
A in Eq. ([SI]) by Ai{a) = 2A/{1 - a). The factor of ^^(0) on the left-hand side of Eq. 
correctly accounts for the normalization of the soft function. (Of course, from Eq. (1641) . we 
see that Sa{0) = 1, but the normalization would not be trivial in the analogous equation 
for the jet function, for instance.) In terms of the matrix elements above, we find a general 
form for the coefficients A„, which is not limited to NLL resummation. 



A„(a) 



1 — a 



An{a) 



(66) 



which, in the general case for n > 1, may still depend upon a, as indicated. 

To explore the information contained in the cumulants. An, let us study the a dependence 
of the parameters A„ in Eq. (!66|) for low n. The n = 1 term, Ai(a), is independent of a, 
as we showed in the previous section, in agreement with the resummed perturbation theory 
result, Eq. (15^ . The a dependence of the second and higher terms, however, differs in 
general. Nevertheless, boost invariance always allows us to perform one rapidity integral in 
the cumulant matrix elements. For the case n = 2, we have 



A2(a) 



dJl + (1 - a) |r7|le-l^l(i-")((fT(0)^T(r7))). 



(67) 
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Under certain conditions, the a dependence of this expression also disappears. In Ref. |16| . 
it was observed that the scahng rule for the nth cumulant term in Eq. (!65|) is good when 
the energy flow correlations are negligible for rapidity intervals larger than a range A?7 ~ 
l/[rz(l — a)]. Assume, then, that the correlator {{St{0)£t{v))) is nonzero only for t] <^ 2{i-a) ■ 
Then we may Taylor expand the remainder of the integrand in Eq. (l67l) about 77 = 0: 



l^dvi^l - - a)vr + ^[(1 - a)v]' + ■ ■ ■ } ((£:t(0)£:H^))). (68) 

Insofar as the correlator {{St{0)St{'>]))) has support only over a region rj <^ the leading 
term of the expansion dominates, and we recover the a- independence of A2. Interestingly, 
there is no 0{ri) term in the expansion multiplying the correlator, so that violations of the 
scaling rule should be even smaller than one might initially expect, at least for moderate 
values of 1 — a. 

We must wait on the analysis of data to interpret the significance of scale breaking for 
the angularities. Supposing, however, that substantial scale breaking were found in the 
power (zz/Q)^ in the shape functions for angularities, we can learn about nonperturbative 
correlations in energy flow through Laplace moments of the cumulants. For example, using 
Eq. ( !67l) . we observe that 

d 

^2(0) - g^^^^ _ ^2(0) = A2(a), (69) 

where C2{a) is a direct Laplace moment of the correlation operators in terms of their rapidity 
separation, with one fixed at rapidity zero, 

C2(a) = - y e-(^-)'^l {{STmrm . (70) 

Assuming that the correlations vanish for a —00, the solution to Eq. fl69l) gives these 
Laplace moments directly in terms of the cumulants A2, and hence in terms of the coefficient 
A2(a), which is, in principle, an observable, 

r da' 

C2(a) = (l-a) / -^A2(a'). (71) 



Furthermore, derivatives of C2{a) with respect to a provide information on Laplace moments 
of the cumulants supplemented by powers. 

For higher n, the situation becomes somewhat more complex, but continues to encode 
potentially interesting physical information. For the coefficient ,43 in Eq. (jHS]), we can 
similarly perform one of the three rapidity integrals and obtain 

Ai{a) = - / rfr/2 / rfr/3e-(^^+''3)(i-a) 

1 - a io io (72) 

X [3((^r(0)^T(-^2)^T(^3))) - {{STmT{r]2)ST{m)))] , 

which again respects the scaling under the assumption that the correlators are nonzero 
only for 772,3 < 
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VII. CONCLUSION 



We have explored power corrections for event shapes using factorization theorems in both 
full QCD and SCET. In this context, we have shown that the formalisms lead to equivalent 
event shape functions that summarize nonperturbative effects of soft gluon emission on event 
shape distributions for two-jet events. We have shown how the boost invariance of lightlike 
Wilson lines implies the universality of the leading 1/Q corrections to the mean values of 
the event shapes, without relying on low-order or even resummed perturbative calculations. 

In addition, we have used the field-theoretic formalism to interpret potential violations 
of the scaling rule for angularity shape functions in terms of correlations between energy 
flow operators for soft gluon radiation. Using l/Q"^ corrections in shape functions as an 
example, we have demonstrated how, in principle, a violation of scaling for the angularities 
can provide information on specific matrix elements in the effective theory. The analysis 
of existing extensive and high-quality data from leptonic annihilation experiments could, in 
this way, provide a new experimental window into the process of hadronization in quantum 
chromodynamics. 
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